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The density matrix renormalization group (DMRG) has been extended to study quantum phase 
transitions on random graphs of fixed connectivity. As a relevant example, we have analysed the 
random Ising model in a transverse field. If the couplings are random, the number of retained states 
remains reasonably low even for large sizes. The resulting quantum spin-glass transition has been 
traced down for a few disorder realizations, through the careful measurement of selected observables: 
spatial correlations, entanglement entropy, energy gap and spin-glass susceptibility, among others. 
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I. INTRODUCTION 

Quantum phase transitions constitute a very active 
topic in theoretical physics^, to which disorder and 
frustration of a quantum spin glass add considerable 
intricacies 2 -. Although quantum effects were usually con- 
sidered to be negligible in spin glasses at finite temper- 
atures, diverse experiments^^ have proved them to be 
highly relevant, leading to the first experimental real- 
ization of quantum annealing. One of the most widely 
employed theoretical approaches to the quantum spin- 
glass transition (QSGT) is the random couplings Ising 
model in a transverse field (RITF) 7 -. Within the ana- 
lytical approach, a Griffiths phase was found in the ID 
case via an elegant RG procedure^, while such a struc- 
ture was proved to be absent from the infinite range case 
within the replica formalis m 10 : 11 : 12 . The ID RG analysis 
was extended to study the presence of a certain num- 
ber of long-distance links, proving the relevance of the 
perturbation^. Regarding numerical approaches, quan- 
tum Monte-Carlo (QMC) remains as the most suitable 
tool at T > 0. Using it, a coherent picture was obtained 
in the 2D and 3D cases, showing that the Griffiths phase 
is present in 2D but absent in 3 D 14 : 15 : 16 . 

Classical spin glasses have been thoroughly analysed 
on Bethe lattices^ 7 -. However, finite Bethe lattices are 
dominated by boundary effects and present no frustra- 
tion. A way out is to study random graphs with a fixed 
connectivity which, in the classical case, allow the appli- 
cation of the cavity approach^. These random graphs 
present genuine frustration, since they have loops, and 
have no boundaries. On the other hand, the average size 
of the loops grows with the number of sitesi^, thus mak- 
ing them resemble, locally, a Bethe lattice (see figure [l}. 

The density matrix renormalization group 
(DMRG) 20 ' 21 is known to be a highly accurate method 
to analyze T = ground state properties of ID or 
quasi- ID quantum many body systems, including tree 
structures^. We have extended the method in order 
to trace the behaviour of finite size samples of random 
graphs of fixed connectivity K = 3 across the QSGT 
within the RITF model. Given the high connectivity 



of these graphs, the applicability of the DMRG is 
highly non-trivial. It is remarkable that in the case of 
non-random couplings, we have found the number of 
retained states needed to ensure a good accuracy to 
be very high, rendering the application of the DMRG 
unfeasible. 

The careful analysis of pseudo-critical points of finite 
size samples constitutes a powerful tool to study a QSGT. 
As a recent example, the work of Igloi and coworkers^ 
studies the distributions of, e.g., the average entangle- 
ment entropy and the surface magnetization of finite size 
ID samples. In this work we provide measurements of, 
among other observables, spin-glass susceptibilities, en- 
ergy gaps, long distance spin-spin correlations and en- 
tanglement entropies, and use them to obtain insight on 
the mechanism of the transition. A full characterization 
of the QSGT, nonetheless, is not attempted in this work, 
since it would require to obtain statistics on a large num- 
ber of disorder realizations. 

The rest of the paper is organized as follows: Sec. |TT] 
presents the model, and Sec. Hill provides details on 
the application of DMRG. Sec. IIVI shows the numerical 
DMRG results obtained. Sec. fVl presents a discussion of 
the results and some concluding remarks. 



II. MODEL 

Let us consider the random Ising model in a transverse 
Field (RITF) on a generic graph 7 -: 
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where denotes pairs of neighboring sites i and j on 
the graph. The values of Jy are uncorrelatcd random 
variables with a uniform probability density distribution 
in the interval [—1,1]. We will focus our study on ran- 
domly generated graphs, with sites and a fixed connec- 
tivity K = 3 18 ' 19 , i.e., each site i is connected to K = 3 
other (randomly chosen) sites j. An example of such a 
graph, with N = 30 sites, is shown in Fig. [TJ 
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We will take the product of eigcnstates of of, for all i, 
as the canonical basis for our problem. In this basis, we 
denote the state in which all spins point in the positive 
x-direction as \x) = 2 _Ar / 2 JTi(IT) + II))- Let us assume, 
for the moment, that h\ = for all i. For T = 0, all the 
non-diagonal terms in the Hamiltonian vanish. There- 
fore, the ground state is given by the configurations with 
the classical minimum energy. Disregarding accidental 
degeneracies, there are two such configurations, related 
by a simultaneous flip of all spins, erf — » —erf for all i. 
The transverse field T may be considered as a kinetic 
energy coefficient, providing a hopping term among the 
classical configurations. In the T — » oo limit the trans- 
verse field term dominates and the ground state tends to 
the state \x), with all the spins pointing in the positive x- 
direction, separated by a large gap A m 2T from the first 
excited state. Let us remark that all the components of 
the canonical basis have the same amplitude in the state 
\x). 

By decreasing T there is a certain value T — T c for 
which the system undergoes a quantum spin-glass tran- 
sition (QSGT). For T > r c , the system is in a quantum 
paramagnetic phase, which presents no long-range order. 
Below r c , the system is in a quantum spin-glass phase, 
presenting a hidden long-range order which may be de- 
tected by a number of observables. We will focus here on 
the divergence of the spin-glass susceptibility: 



FIG. 1: (a) A random graph with N — 30 sites and connec- 
tivity K = 3. Stronger bonds highlight the local Bethe lattice 
structure around site 25. Numerical labels denote the order in 
which they are visited along the DMRG path, (b) Deforma- 
tion of the same graph to show the DMRG path as a straight 
line. The left and right blocks at a certain DMRG step are 
shown, and the dashed line cuts all the long-range links that 
would need reconstruction during that DMRG step, 9 in the 
example. 



It should be noticed that, locally and for large TV, such 
graphs resemble Bethe lattices, as we have highlighted 
with thicker links in Fig. Q] Nonetheless, the behaviour 
of finite Bethe lattices is dominated by the boundary, un- 
like in our case, where there is none. Moreover, trees have 
no loops and, therefore, present no frustration, while ran- 
dom graphs with fixed connectivity do contain loops and, 
therefore, have genuine frustration. It has been proved, 
nonetheless, that short loops become rare as TV — » oo: 
More precisely, the probability of finding a loop of any 
fixed size L falls to zero when TV — > oo^. The physics 
of classical spin glasses in these graphs has been studied 
using the cavity method by Mczard and Parish. 

Let cat be the absolute value of the disorder-averaged 
energy per site of the classical ground state. Our nu- 
merical experiments show that it increases from ~ 0.6 at 
TV = 100 up to ~ 3/4 for larger TV, which is the theoreti- 
cal limiting value for TV — > oo under the assumption that 
all links arc satisfied. 



XSG = -77 



1 V /d(crf) 
TV ^ l dh z A 



(2) 



i.e., physically, a small longitudinal magnetic field h*-, 
applied at site j, generates a magnetization response on 
each site i, which is measured (and squared, so as to dis- 
regard its sign); the results arc summed over all sites i 
and averaged over all sites j. If the system is in a para- 
magnetic phase, the magnetization will be proportional 
to h* and short-ranged in space, so that the sum over i 
and j yields a finite value for xsg- On the other hand, 
on approaching the spin-glass phase, an infinitcsimally 
small longitudinal magnetic field /i| , localized at a single 
site j, will eventually induce a finite magnetization over 
a long-range of spins. This effect is at the origin of the 
divergence of XSG- 

We shall now discuss the numerical approach we have 
used to study this system, and the results obtained. 



III. APPLICATION OF THE DMRG 

The density-matrix renormalization group (DMRG) 
has proved to be an accurate method to analyse the prop- 
erties of ID and quasi-lD systems^-^. The DMRG 
may be described as a variational method within the 
matrix-product states (MPS), which constitute a low- 
dimensional subspace of the full Hilbert spaced. A MPS 



3 



may be expressed as 

|*>= J2 Tr (a (3i) • • • A^) \ Sl ---s N ) , (3) 

si---sjv 

where each A^ Si ' is a square matrix with dimension m, 
which may be considered as the number of retained states 
per block when splitting the system into a left and right 
parts. The total number of variational parameters is less 
than 2Nm 2 < 2 N . The success of the DMRG is related 
to the ability of these MPS to reproduce faithfully the 
ground states of local ID many-body problems for low 
values of nv&. If m — ► oo, any state of the Hilbert space 
may be exactly represented as a MPS. 

The minimum number of retained states m is related 
to the exponential of the von Neumann entanglement en- 
tropy of the DMRG blocks. In a non-critical ID system, 
this entropy is bounded for all sizes, while it grows as InL 
for a ID critical system of length L. It is believed that 
for a D-dimensional system out of criticality, the entan- 
glement entropy scales as L , where L is the short- 
est spatial dimension of the system^. This estimate is 
known as the area law and is believed to have logarith- 
mic corrections at critical points. An important practical 
consequence is that, in order to study a 2D system, the 
number of retained states m should grow as exp(L), thus 
making the DMRG very inefficient. 

Our system, on the other hand, is defined on a ran- 
dom graph of fixed connectivity, K = 3. We will show 
that this poses no problem to the number of retained 
states m, which appears to remain manageable even for 
N = 500 as long as the couplings are random. However, 
the implementation of the DMRG on such a model has 
required the following technical refinements of the origi- 
nal method: 

(a) Path selection. In a non-ID system, DMRG pro- 

ceeds by converting the system into an effective ID 
problem with long-range couplings. A path is cho- 
sen along the graph, which does not repeat sites, 
and is considered to be appropriate if the num- 
ber of broken links between the left and the right 
blocks is kept low along a DMRG sweep. Normally, 
the selection of a suitable path in a quasi-lD sys- 
tem (e.g. ladders) is done by geometrical intuition. 
In our implementation we have designed an auto- 
mated procedure: a simulated annealing algorithm 
is employed in order to minimize the number of 
broken links. The full problem of finding the opti- 
mal path is computationally very hard. Therefore, 
we do not aim at the exact optimum, but only to a 
reasonably good local minimum. We have noticed 
that our quasi-optimal path performs much better 
than a random path. 

(b) Perron-Probenius criterion. The Hamiltonian of 

the RITF on any graph fulfills, on the canonical 
basis, the conditions of the Perron-Frobenius the- 
orem, i.e., all off-diagonal components are non- 
positive. Hence, all the ground state components 



must have the same sign. Using the MPS repre- 
sentation of the ground state obtained within the 
DMRG, it is always possible to reconstruct the am- 
plitude of any configuration C = {s\ ■ ■ ■ sn}. The 
obtention of the full \^>) is unfeasible, since it would 
require reconstructing the amplitude of an expo- 
nentially large number of configurations. Nonethe- 
less, it is possible to pick up a few random configu- 
rations C, and check that all their amplitudes have 
the same sign. Whenever this criterion was not 
met — a rather rare event — , the calculation was 
repeated changing the random seed for the Lanc- 
zos procedure on the first DMRG step. 

(c) Wavefunction annealing. For large T the ground 

state \x) is easily represcntablc as a MPS, requir- 
ing a single retained state. The DMRG works ex- 
tremely well in this regime and, therefore, our sim- 
ulations are always started well within the param- 
agnetic (large T) phase. The QSGT is approached 
by repeatedly decreasing the value of T by a small 
amount, always using the previous ground state as 
a seed for the new calculation, exploiting the wave- 
function transformations suggested by White^. 
This type of annealing of the ground state provides 
a faster convergence and more accurate results for 
low r. 

(d) Adaptive number of retained states. The 

number of retained states m, and the number of 
DMRG sweeps n s , are not fixed in our algorithm. 
We set a maximum value for the sum of the 
neglected eigenvalues of the density matrix in the 
RG truncation « 10~ 6 ), and select m accord- 
ingly. Moreover, sometimes convergence takes 
more sweeps than usual (n s ~ 30 — 40) in order to 
obtain machine precision in the convergence of the 
ground state energy. 

(e) Energy gap measurements. When there is a 

symmetry in a problem, e.g., under SU(2), it is 
usually possible to obtain the first excited state as 
the ground state of a different sector of the Hilbert 
space. In our case, lacking this, the best option has 
proved to be the following one. For each DMRG 
step, after the ground state |^o) had been found, 
it was "promoted" to a higher energy by the fol- 
lowing transformation of the Hamiltonian: 



H H + A|* > (*o| , (4) 

in such a way that the ground state of this new 
hamiltonian is the former first excited state, as long 
as we take A > A, A being the gap we are looking 
for. The density matrix used for truncation was 
built as a linear combination of the density matri- 
ces for each state, with equal weights. However, we 
should remember that only ground states of local 
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Hamiltonians are expected to be faithfully repre- 
sented by MPS2£, and Eq. j4} does not define a 
local Hamiltonian. Therefore, the accuracy in the 
gap estimate is worse than that in the ground state 
energy and in other observables. 

IV. RESULTS 

We have applied the modified DMRG technique to 
study a few samples of random graphs with fixed connec- 
tivity K = 3 and sizes ranging from N = 50 to TV = 500. 

A full characterization of the QSGT would require rel- 
evant statistics on the disorder. Unfortunately, in order 
to obtain a high accuracy for each sample, the required 
CPU-time is rather large. Therefore, we have decided to 
focus on a precise study of a few realizations, in order 
to gain insight on the mechanism of the transition for 
finite samples. Most of our findings will be illustrated in 
the figures of this section by showing in detail the results 
obtained for a sample with N — 200 sites. It should be 
remarked that the highlighted features are typical of all 
the ensemble. 

The numerical simulations were done with the DMRG 
algorithm explained in section HTH with a neglected prob- 
ability tolerance r\ = 10~ 6 and a tolerance on the conver- 
gence of the energy of one part in 10 10 , for ten samples 
of each size (N = 50, 100, 150, 200, 250, 300, 400, 500). 
Each sample is a different random graph, always with 
fixed connectivity K = 3, and the bonds are random 
and independent, uniformly distributed in the interval 

[-1,1]- 

Figure [2ja) shows three observables calculated for a 
given instance of a N — 200 graph, as a fuction of T. The 
most relevant quantity is the spin-glass susceptibility, de- 
fined in equation [2 A small magnetic field h\ = 10 -4 is 
applied at a single site io, and the magnetic response is 
measured with the formula 

which is found to increase very fast as T approaches 
1.28 from above, and thereupon saturating at a very 
high value. Figure HJb) proves the divergence of XSG 
by showing its behavior at three different values of h\ 
(h* = 10~ 3 , 10~ 4 and 10~ 5 ): the saturation value is 
seen to scale as ) 2 . We estimate the pseudo-critical 

value r c for the given sample as the value of T at which 
the slope of the susceptibility attains its maximum. With 
this definition, r c appears to be almost independent of 
the chosen site io, at our level of precision AT = ±0.01. 
Figure [2^c) , finally, shows xsg for ten different samples 
with N = 200 sites, differing in the graph structure and 
in the choice of the couplings Jy , and diverging at dif- 
ferent (sample dependent) values of T c . 

Figure [Ha) also shows two other observables which 
point towards the same value for T c . The first is the block 



entanglement entropy, defined as S = — Tr(plog(p)), 
where p is the reduced density matrix for a part of the 
system. It is measured for all the left-right decomposi- 
tions along a DMRG sweep, and its maximum value is 
denoted by S max ■ This value of S mnl is obviously depen- 
dent on the DMRG path. Nonetheless, since this path 
has been chosen so as to minimize the number of retained 
states, it is expected that it will also minimize the maxi- 
mum value of the entropy. It should be emphasized that 
this observable attains its maximum value, as a function 
of r, at the same T c which is found by analyzing the 
spin-glass susceptibility. The relevance of the entangle- 
ment entropy in order to characterize the critical point of 
a QSGT has already been remarked in the literature^. 
The DMRG, being based on the obtention and analysis 
of the density matrix of different parts of the system, is 
specially well suited for measuring this observable 2 ^. 

The other observable shown in figure EJa) is the energy 
gap A, which extrapolates to zero at a value of T which 
is close to r c . The sign of the difference between these 
two pseudo-critical points is sample dependent. This dis- 
crepancy is likely to be a finite size effect. 

The next observable that we calculated for each sample 
is the average value of (S x ), sec figure [3] For high T, the 
transverse field is dominant, the ground state is close to 
\x) and the value of (S x ) is close to 1. As we decrease T, 
this value is reduced. It is remarkable that, at the value 
of r c defined by xsg> an d confirmed by the maximum 
block entropy, (S x ) is never lower than 0.95, at variance 
with the ID RITF, where it is normal to have values be- 
low 0.5. Also in figure [3] we have plotted the behavior 
of the quantum equivalent of the q^ order parameter, 
which measures the magnitude of the long-distance spa- 
tial correlations^, 

Our definition differs slightly from that used commonly 
in the literature on classical spin glasses. Instead of the 
long-distance behavior, we measure the global average 
behavior. They are equivalent in a system with a Bethe 
lattice topology because the number of neighbours of a 
given site at distance d scales as exp(d). Also note that 
these two observables, (S x ) and q^ have the advantage 
that they do not vanish trivially in the absence of external 
longitudinal magnetic field. 

This quantity is close to zero within the paramagnetic 
phase, and takes a non-zero value inside the spin-glass 
phase. Figure [3] shows that the behaviours of q^ 2 ' and 
(S x ) are highly correlated, both pointing to a pseudo- 
critical point which is, in many instances, sensibly lower 
than the one obtained with XSG and S max . This sys- 
tematic discrepancy is not explainable in the framework 
of classical physics, since the fluctuation-dissipation the- 
orem leads to the proportionality of g' 2 ' and xsg^- In 
the quantum case, XSG is not equivalent to q^ , but to 
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FIG. 2: (a) For a random sample with N — 200 sites, the 
spin-glass susceptibility xsg is shown in logarithmic scale, 
along with the maximum block entropy 5 m „ and the energy 
gap A, the two latest in arbitrary units. The divergence of 
Xsg is marked with the vertical line, which clearly coincides 
with the maximum value of S m ax ■ The energy gap comes very 
close to zero (dotted line) precisely in that region, (b) The 
fact that xsg saturates is related to the finite value of the 
applied h z field. In this figure we can see how this saturation 
value increases as h z is decreased, (c) The divergence of xsg 
is shown for several samples, all with N = 200 sites. 
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FIG. 3: The behavior of (S%) and the g*- 2 ' order parameter is 
shown for the same sample of figure [2] The two curves are 
strongly correlated, suggesting a transition at the same value 
of r. This value is substantially lower than the F c pointed by 
Xsg, and S m ax, which is marked in the plot with a vertical 
line. 
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(7) 



where the integral is performed on imaginary time, from 
t = to j3. Therefore, in quantum spin-glasses, xsg 
contains a contribution from time-correlations, while 
only measures spatial ones. Their different divergence 
points might suggest that long-distance correlations de- 
velop at a higher value of F in time than in space. 

At different sizes TV, the behavior of the various sam- 
ples is qualitatively similar, the only difference being the 
sample-dependent value of T C (N). It is possible to con- 
struct histograms showing, for different N, the probabil- 
ity distribution of the various T C (N). As a preliminary 
result, we mark in figure [4] by empty circles the r c values 
obtained for 10 samples of different sizes: A^=50, 100, 
150, 200, 250, 300, 400, and 500. The crosses mark, for 
each size N, the average values of T C (N), which appar- 
ently saturate at some point around 1.5. 

A useful clue to the physics of this system across the 
transition is obtained by monitoring selected wavefunc- 
tion components of the ground state as a function of T, 
which can be easily done with the DMRG. We illustrate 
this in figure El for a sample with TV = 100 whose T c is 
marked by the arrow. One of the monitored configura- 
tions is the classical ground state of the system. (Natu- 
rally, there are two of them, related by a global spin-flip, 
which we denote by <f> + and (f>~ . They are obtained by 
measuring the values of (Sj) at each site for very low 
r under the presence of a very small longitudinal field 
h z that splits the degeneracy between them. We will 
consider its weight to be the sum of their probabilities.) 
The weight of the classical ground state grows up to 1 



FIG. 4: The abscissa of each circle denotes the value of r c for 
a different sample, its ordinate being its size N. The crosses, 
joined by a line, denote the average values of T C (N) for each 
size. These average values seem to converge to a value around 
1.5. 

as r is decreased. A second monitored configuration is 
picked at random (dashed lower line in figure [5]) : the 
weight of such a random configuration is similar to that 
of all the others only for very large values T, while it is 
definitely much smaller when T approaches the pseudo- 
critical point. The remaining monitored configurations 
in figure [5] are obtained by classical simulated annealing, 
i.e., they are local minima of the classical energy. These 
configurations maintain a high weight (similar to that of 
the optimal configuration) across the transition, up to a 
value of r below which their weight decline markedly. 

Therefore, the number of relevant classical configura- 
tions contributing to the ground state changes drastically 
across r c . Deep into the paramagnetic phase, all of them 
have the same weight, coherently bound within the \x) 
state. This state retains a high weight at r c , as shown by 
the high values taken by (S x ) at that moment. Within 
the spin-glass phase, the weight of the different configu- 
rations is gradually redistributed according to their clas- 
sical energies, until eventually only the classical ground 
state remains. 

It should be remarked that the DMRG was applicable 
in practice only in the case of random couplings. For 
random graphs with fixed values of J (cither fcrro or 
antifcrromagnetic), the number of retained states needed 
for a similar accuracy increased in an order of magnitude, 
rendering impractical the calculations. 



V. DISCUSSION AND CONCLUSIONS 

It should be remarked that a quantum phase transition 
may be observed exactly in the MPS formalism with as 
few as two retained states per block^i. The low values 
of S max measured in our system points to the fact that 
the number of states involved in the QSGT is reduced. 



FIG. 5: Probabilities of various classical configurations within 
the ground state, as a function of Y, for a sample with N = 
100 sites. Notice the log scale. The QSGT as obtained from 
the divergence of the spin-glass susceptibility is marked with 
an arrow. One of configurations is the one which minimizes 
globally the classical energy (for F = 0). Another one is 
chosen at random. The rest are configurations with energies 
very close to that of the global minimum. It is noticeable how, 
after the transition, all the states but the random one increase 
their probabilities. Nonetheless, after some further decrease 
in r, all of them but the real minimum reach a maximum 
value and eventually fall to zero. 

Thus, in an attempt to explain the transition we might 
build a simple-minded Ansatz consistent of a linear com- 
bination of only two states: the classical ground state 
\<j>) = (\4> + ) + \<p~)) and a background state containing 
the rest of classical configurations: 

lB) = 7W=2 £ l*>- ^ 

Neglecting matrix elements which are exponentially 
small for large N, there is a very sharp crossover at 
r c (iV) = €n, from state \B) to state |^). This last state 
clearly presents a divergent spin-glass susceptibility. The 
physics of our model is actually more complicated than 
this simple Ansatz, as shown by the numerical estimate, 
T C (N) ps 2e w for large TV. 

An interesting question is how our system differs from 
the ID RITF analyzed by Fisher and coworkers^. In 
ID, duality arguments and a detailed RG analysis give a 
value of T c satisfying logr c = [log Jij] a v, where [-] av indi- 
cates a disorder average. This equation leads to values of 
r c much lower than those measured in our random graph 
case (e.g., for our distribution of Jy, Fisher's model 
yields T c = e _1 , while in our case it reaches r c ps 1.5 
for N = 500 sites). The higher connectivity seems to 
make a large difference in that respect. Moreover, in our 
case, the value of (S x ) at the transition is fairly high, 
about 0.95, while it is almost always below 0.5 for RITF 
chains with the same sizes. This means that the state 
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\x) is still dominant in the paramagnetic phase at the 
moment of the divergence of the spin-glass susceptibility. 

In conclusion, we have extended the DMRG to make 
it suitable to study the QSGT on random graphs. The 
main technical innovations are the path-selection, which 
reduces the maximum number of retained states, and the 
wavefunction annealing, which allows to use the ground 
state for a certain value of T as a seed to obtain the 
ground state for a lower value. 

This modified DMRG algorithm has been applied to 
the measurement of the energy gap A, spin-glass suscep- 
tibility xsGi entanglement entropy S max , long-distance 
spatial correlations and average x-magnetization 
(S x ) to high precision on a few samples of different sizes, 
ranging from N = 50 to N = 500 sites. Remarkably, 
Smax attains its maximum at the same value of T = T c 
at which XSG diverges. This led us to consider this value 
as our candidate for the pseudo-critical point. The en- 
ergy gap A vanishes in the surroundings of that value for 
all realizations. In most disorder realizations, and 
(S x ) start their crossover at a value of T which is inferior 
to T c . This may suggest that long-range temporal corre- 
lations develop at a higher value of T than purely spatial 
ones. 

We would like to emphasize that our definition of S m ax 
constitutes an alternative approach to the entanglement 
entropy on random graphs. Its relation to the standard 
block entropy, defined as the entanglement entropy ob- 
tained by tracing out a site and its neighbours up to a 



certain distance, remains as an open question. 

We have also monitored the ground state probabili- 
ties of several classical low-energy configurations, show- 
ing that, as r is reduced below r c , these probabilities in- 
crease exponentially until they attain a maximum value 
and then fall to zero, leaving the classical minimum en- 
ergy configuration as the only ground state component 
as r -> 0. 

The fact that the model on such disordered graph was 
amenable to analysis within the DMRG is non-trivial. 
Both S m ax and the maximum number of retained states 
per block increase slowly with the system size. In a re- 
stricted sense, the system behaves similarly to a ID chain, 
despite the high connectivity of the underlying graph. 
This is a rather remarkable effect due to the disorder, 
which perhaps leads to the selection of an effective ID 
path of strong bonds. If the couplings J are not random, 
we have found that the number of retained states per 
block increases in a much more pronounced way, render- 
ing the numerical DMRG approach impractical. It will 
be interesting — and is left to a future study — to under- 
stand quantitatively how the entanglement entropy at 
the transition behaves as a function of the system size, 
both for random and non-random coupling o 32 ' 33 ' 34 . 
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